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Abstract — We clarify from a general perspective, 
the condition for the appearance of chaotic energy 
spectrum in quantum pseudointegrable billiards with 
a point scatterer inside. 

I. Introduction 

The quantum billiard with point scatterers is a quasi- 
exactly solvable model which is closely related to real 
systems. The billiards are a natural idealization of 
the particle motion in bounded systems. The single- 
electron problem in mesoscopic structures is now a 
possible setting, owing to the rapid progress in the 
mesoscopic technology. Real systems are, however, 
not free from impurities which affect the particle mo- 
tion inside. In the presence of a small amount of con- 
tamination, even a single-electron problem becomes 
unmanageable. The modeling of the impurities with 
point scatterers makes the problem easy to handle 
without changing essential dynamics. 

A fundamental problem for the billiards considered 
here is to understand global behavior of the energy 
spectrum in the parameter space of particle energy 
and the strength of the scatterers. In particular, sta- 
tistical properties of the spectrum are important be- 
cause they reflect the degree of complexity of underly- 
ing dynamics. It is widely believed that integrable sys- 
tems obeys Poisson statistics, while the predictions of 
the Gaussian Orthogonal Ensembles (GOE) describe 
chaotic systems (I]] . 

In this paper, we discuss the spectral properties of 
quantum billiards with a single point scatterer. It 
should be noticed that the nature of classical motion in 
such systems is pseudointegrable Q in the sense that 
unstable trajectories, which hit the point scatterer, are 
of measure zero. As a result, it is expected that the 
energy spectrum of a quantum analogue does not sub- 
stantially differ from Poisson statistics. However, as 
shown in this paper, quantization induces the chaotic 
spectrum under a certain condition. This phenomenon 
might be called wave chaos || because its origin is the 
wavelike nature of quantum dynamics. After deriving 
the eigenvalue equation with the help of self-adjoint 
extension theory in functional analysis, we discuss the 
statistical properties of the energy spectrum from a 
general perspective and deduce the condition for the 
appearance of chaotic spectrum in the quantum pseu- 



dointegrable system. 

II. Formulation 

We start from an empty billiard in spatial dimension 
d, d — 1,2,3. Let us consider a quantum point parti- 
cle of mass M moving freely in a bounded region 
We impose the Dirichlet boundary condition such that 
wave functions vanish on the boundary of f!^. The 
eigenvalues and the corresponding normalized eigen- 
functions are denoted by E n and ip n (%) respectively; 



(1) 



The Hamiltonian Ho is the kinetic operator in L 2 (fi^^) 
with domain D(H ) = H 2 (ttW) n Hb(nW) in terms 
of the Sobolev spaces. The Green's function of Hq is 
given by 



G(°)(£,f;a,) = £ 



n=l 



lf n (x)ip„(x') 
~ En 



(2) 



The average level density is given by 
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d= 1, 
d = 2, 
d = 3, 



(3) 



where we denote the "volume" of fl^ by the same 
symbol. (For example, the area of a two-dimensional 
region is simply denoted by Q( 2 \) 

We now place a single point scatterer at xo inside 
the billiard. Naively, one defines the scatterer in terms 
of the d-dimensional Dirac's delta function: 



H = H a + vS (d) (x~ x ). 



(4) 



However, the Hamiltonian H is not mathematically 
sound for d > 2. It is easy to see that the eigenvalue 
equation of H is reduced to 



E 



LU - E n 



(5) 



Keeping Eq.(H) in mind, we realize that the infinite 
series in Eq.(g) does not converge for d > 2, since the 
average of the numerator (among many n) is energy- 
independent; (fn(xo) 2 ) — 

To handle the divergence, a scheme for renormaliza- 
tion is called for. One of the most satisfying schemes 
is given by the self-adjoint extension theory of func- 
tional analysis jjj. We first consider in L 2 {Q,^) the 
nonnegative operator 



where 



2M 



(6) 



with its closure Hg in L 2 (Q^). Namely, we restrict 
D(Hq) to the functions which vanish at the location 
of the point scatterer. By using integration by parts, 
it is easy to prove that the operator Hg is symmetric 
(Hermitian). But it is not self-adjoint. Indeed, the 
equation 

Hi^ u [x)=^ u [x), ^eD(Hl) (7) 
has a solution for Im u> ^ Q ; 

iP u (x) = GW(x,x ;u), xen^-xo (8) 
indicating 

D(Ht Q ) = D(H So ) © N(H* So - w) © N(Hl - Co) 

^D{H So ), (9) 

where N(A) is the kernel of an operator A. Since 
Hg has the deficiency indices (1,1), Hg has one- 
parameter family of self-adjoint extensions Hg ||; 

D(Hg) = {/|/ = <p + c(iP iA - e 4 V-^A); 

^D(4),ceC,0<ff<2f}, 
Hgf = H So cp + iAc(ip iA + e l V_ 4A ), (10) 

where A > is an arbitrary mass scale. With the 
aid of Krein's formula, we can write down the Green's 
function for the Hamiltonian Hg as 

Gg(x,x,Lo) = G w (x,x';Lo) 

+ G^(x,x Q ;u J )Tg(io)G (a \x ,x';u J ). (11) 

The transition matrix (T-matrix) Tg is calculated by 



T e (u) 



where 



jo 



(lo — iK)ci A {uj) — e lS (uj + iA)c-iA(u) 



(12) 



c ±iA (u;) = / G^{x 1 x ;Lo)G^{x,x ;±iA)dx. (13) 
Jaw 



Using the resolvent equation, we have 
T g (u) = (vZ 1 - G(w)) _1 , 



(14) 



-i = A cot^y^#^, 



oc 1 

GM = ^^„(f ) 2 ( 



Ey, 



n=l 



,J - E n El + A 2 



(15) 



)• (16) 



The constant vg is a coupling constant of the point 
scatterer, the value of which ranges over the whole 
real number as one varies < 8 < 2tt. It follows from 
Eq.(fi3) that the eigenvalues of Hg are determined by 



G(co) 



(17) 



On any interval (E n , E n +i), G is a monotonically de- 
creasing function of lo that ranges over the whole real 
number. This means that Eq.(J17j) has a single solution 
on each interval for any vg. The eigenfunction of Hg 
corresponding to an eigenvalue u> n is given by 



{X,X ;U) n ). 



(18) 



III. Condition for Strong Coupling 

A. d = 1 

In one dimension, each infinite scries on RHS in 
Eq.(|l^) converges separately. Thus, if we set 



n=l 



Wn(xoY 



En 



El + A 2 ' 



(19) 



Eq. ( |rl| ) is reduced to Eq. (||) ; A somewhat complicated 
argument above leads us to the well-known result with 
Dirac's delta. In the following, we consider Eq.(|^), the 
LHS of which is denoted by G (1) M- 

To obtain each solution of Eq. (|B|) , a numerical task 
is needed in general. However, since our main pur- 
pose lies in the statistical properties of spectrum, we 
proceed further by introducing some approximations 
without losing the essence, while still keeping loss of 
generality minimal 0. The first is that the value of 
<fn(xo) 2 is replaced by its average among many n; ^ 



Vn {x ) 2 ~ (^„(f ) 2 > = l/« (d) - 



(20) 



Since the statistics is taken within a large number of, 
sometimes thousands of eigenstates, Eq.(|o[) is quite 

satisfactory. Keeping [(gnfafn)] — l/vW) in mind, we 
recognize from Eq.(|18|) with Eq.(|2|) that if u> m ~ E m 
(or E m+L ) for some to, then ip m ~ (p m (or ip m +i), 
implying that only ip m with an eigenvalue u> m ~ (E m + 
E m+ i)/2 is distorted by a point scatterer. For such 
to m , since the contributions on the summation of G^ 
from the terms with n ~ to cancel each other, G^ 
can be estimated by a principal integral; 



En 



E n 



(21) 



1 With some minor modifications if necessary, a major part 
of the argument here is applicable to higher dimensions. We 
thus explicitly signify u d = 1" only for the assertions specific to 
one dimension. 



o w-E 



dE, d=l, (22) 



indicating that the summation in Eq.(p9|) converges 
rapidly. From Eqs.(||) and ( p0| ) for d = 1, we obtain 



where we have defined a continuous function g( d > of 
to which behaves like an interpolation of the inflection 
points of G (d) . Inserting Eq.|§ into Eq.(||) and using 
the elementary indefinite integral 



a (1) K0 



fM 1 / 2 



2V2 



(30) 



dE 



An 



(23) 



(uj - e)Ve 

for uj > 0, we obtain for u) m ~ (E m + E m+ i)/2 

G {d) (uj m ) ~0, d=l, (24) B. d = 2, 3 



which is inversely proportional to square root of the 
energy uj. We can summarize the finding in one di- 
mension as follows; The effect of a point scatterer of 
coupling strength v is substantial only in the eigcn- 
states with eigenvalue uj such that 



A«(a 



ttM 1 / 2 1 
~ 1 2?^7fu\ 



(31) 



indicating that the maximal coupling of a point scat- 
terer is attained when the strength v satisfies 



0. 



1, 



(25) 



which is going well with our intuition. 

The "width" of the strong coupling region (allowable 
error of u _1 in Eq.(p5|)) is estimated by considering 
a linearized eigenvalue equation. Expanding G^' at 
ui m ~ (E m + E m+ i)/2, we can rewrite the eigenvalue 
equation as 



G^'(uj m )(uj - uj m ) ~ v- 1 - G {d) {uj m ). 



(26) 



In order that Eq.(|2q) has a solution uj ~ uj m , the range 
of RHS has to be restricted to 



v- x -G^^ m ) 



< 



A^(uj m ) 



AW(o 



(27) 



(28) 



where we have defined the width A^ d ' which is nothing 
but the average variance of the linearized G^ d ' on the 
interval (E mi E m +i). Using the approximation (|2o|), 

the value of G^' (uj m )\ can be estimated as follows; 

7^1 V^m - E n 

oo 

- (^n(^o) 2 )^ 



G (<0 '(«m) 



= ^ n (x,f)p^{uo m fY^ 

n— 1 ^ ' 

= n 2 ( Vn (x,f)p[ d }(uj m ) 2 . (29) 

The second equality follows from the approximation 
that the unperturbed eigenvalues are distributed with 

a mean interval pit^(wm)^ 1 in the whole energy re- 
gion. This assumption is quite satisfactory since the 
denominator of G^ (uj) is of the order of (uj — E n ) 2 , 



It is easy to apply the argument above to higher di- 
mensions. Since each series on RHS in Eq.([l6]) diverges 
when summed separately for d — 2,3, we have to re- 
sort to Eq.(p/7J). As a result, Eq.(E2) is replaced by 



<) 



WHH^„(x ) 2 )p 



1 



E 



uj-E E 2 +A 2 
xp^(E)dE, d = 2,3. 
Using the elementary indefinite integrals, 

1 E \ .„ . \u)-E\ 



E E 2 +A- 



E 



dE 



In 



y/E 2 + A 2 ' 



(32) 



(33) 



E E 2 + A 2 



EdE = In 



V^ + Ve 



1 A f E + \/2AE + A 
2V 2 n I £- V2AE + A 



y / ^ jarctan( V ? ^f -!- J ) arctanl \ / I.) 



for uj > 0, 
we obtain for ui m — (E m + E m +i)/2 

d = 2 



(34) 



M 



G"»u;.,)^\ 2 V/"{v^ d = 3 



(35) 



The width is estimated by Eq.(g8|) with Eq. 
Using Eq.(||), we find 




d = 2, 
Wim d = 3. 



(36) 



From Eq.(27) with v replaced by vg, we realize that 
the strong coupling is attained under the condition 



Vn + 



M , uj 

In -r 

2tt A 

M 3/2 A l/2 



2tt 



< 



< 



~4~' 

2 3 / 2 



d = 2, (37) 
VoJ, d = 3. (38) 



IV. Numerical Example 

We consider a quantum particle with mass M = 1/2 
moving in a three-dimensional rectangular box with 
side-lengths (l x ,l y ,l z ) = (1.047,1.186,0.8049) and 
hence il^ = 1. The mass scale is set to A = 1. 
A single point scatterer is placed at the center of 
the billiard. We take into account only the states 
with even parity in each direction, since the others 
are unaffected by the scatterer in this case. Fig.l 
shows, for Vg 1 — 0, 10 and 30, the nearest- neighbor 
level spacing distribution P(S), which is defined such 
that P(S)dS is the probability to find the spacing be- 
tween any two neighboring eigenstates in the inter- 
val (S, S + dS) . (The average spacing is normalized 
to one.) Generic integrable systems obey the Pois- 
son distribution P(S) = cxp(— S) (dotted line), while 
chaotic systems are described by the GOE prediction 
P(S) = irS/2 x exp(-7rS 2 / 4 ) (solid line). The statis- 
tics is taken within umn ~ W3100 in Fig.l (21600 = 
8304). According to Eq.(38), the most strong coupling 



is attained for v e L = -M 3 / 2 A 1 / 2 /2n = -0.05626. As 
Vg 1 increases, P(S) tends to approach the Poisson pre- 
diction. For v^ 1 ~ A( 3 Vi 600 )/2 = 11.3, P{S) is ex- 
pected to be intermediate in shape between Poisson 
and GOE. The value can be considered as the upper 
bound of Vg 1 for inducing a GOE-like shape in this en- 
ergy region. With Vg 1 beyond the bound, the system 
is not substantially different from the empty billiard, 
and as a result, P(S) resembles to the Poisson distri- 
bution. These features are successfully reproduced in 
Fig.l. 

V. Conclusion 

We have discussed the condition for the appearance 
of chaotic spectrum in the quantum pseudointegrable 
billiards with a single point scatterer inside. Chaotic 
spectrum appears if the conditions j3l|), fl37| ) and ( |38| ) 
are satisfied for dimension d — 1,2,3, respectively. In 
two dimension, the condition is described by a log- 
arithmically energy-dependent strip with a constant 
width, indicating that the system recovers integrabil- 
ity in the high energy limit for any vg. In three di- 
mension, the spectrum shows GOE-like nature when 
v^ 1 is within a band whose width increases paraboli- 
cally as a function of the energy. This implies that the 
spectrum becomes chaotic at high energy for any vg, 
which makes a clear contrast with two dimension. 
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Fig.l. P(S) for Vg 1 = 0,10,30. The solid (dotted) 
line is GOE prediction (Poisson distribution). 



